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We report a fully inclusive measurement of the flavor changing neutral current decay 6 — > 37 in 
the energy range 1.8 GeV < E* < 2.8 GeV, covering 95% of the total spectrum. Using 140 fb~^, 
we obtain B{h sj) = (3.55 ± 0.32 1 q - 07) ^ 10"*, where the errors are statistical, systematic 
and from theory corrections. We also measure the first and second moments of the photon energy 
spectrum above 1.8 GeV and obtain (E-y) = 2.292 ± 0.026 ± 0.034 GeV and (E^) - {E-,f = 0.0305 ± 
0.0074 ± 0.0063 GeV^ where the errors are statistical and systematic. 

PACS numbers: 13.20.He, 13.40.Hq, 14.40.Nd, 14.65. Fy 



The flavor changing neutral current decay 6 — s- S7 
process is of remarkable theoretical interest. Its total 
branching fraction is very sensitive to physics beyond the 
Standard Model as it may be affected by the presence of 
charged Higgs or SUSY particles in the loop. Yet the 
present theoretical prediction for the branching fraction 
of (3.791 ^j^) X 10-* jlj,!^, and the average experimen- 
tal value (3.3 ± 0.4) x IQ-^ agree well. This agree- 
ment sets a strong constraint on, e.g., models Q that 
accommodate the observed difference in CP asymmetries 
in the B J/il^Ks and B 4>Ks decays 0. To ob- 
tain stronger constraints on physics beyond the Standard 
Model, more precise theoretical predictions and experi- 
mental measurements are needed. 

On the other hand, the photon energy spectrum 
is almost insensitive to physics beyond the Standard 
Model 01 . At the parton level, the photon is monochro- 
matic with energy E « mb/2 in the b quark rest frame. 
The energy is smeared by the motion of the b quark in- 
side the B meson and gluon emission. A measurement 



of the moments of this spectrum allows a determination 
of the 6-quark mass and of its motion. This information 
can then be used to extract the CKM matrix elements 
\Vcb\ and \Vub \ from inclusive semileptonic b decays 13,111. 
However, a measurement of the low-energy tail of the 
photon spectrum is important in this context 0. 

Belle has previously measured the 6 — > 57 branch- 
ing fraction with 5.8 fb""'^ of data using a semi-inclusive 
approach. Because we applied an effective cut E-y > 
2.24 GeV in the B rest frame, the precision of that mea- 
surement is limited by theoretical errors due to the ex- 
trapolation to the whole energy spectrum. More recently 
the CLEO collaboration has reported a measurement of 
the branching fraction and the energy spectrum moments 
performed in a fully inclusive way [lOj for the range 
E* > 2.0 GeV in the center-of-mass frame Here 
we present a measurement using a similar approach, but 
based on a much larger dataset allowing a detailed study 
of the backgrounds. We extend the photon energy range 
to E* > 1.8 GeV, covering almost the entire spectrum. 
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The b ^ decay is studied using the Belle detector at 
the KEKB asymmetric e^e~ storage ring The data 
consists of a sample of 140 fb^^ taken at the T(4S') reso- 
nance corresponding to (152.0 1 o?) ^ 10^ BB pairs. An- 
other 15 fb~^ sample has been taken at an energy 60 MeV 
below the resonance and is used to measure the non-BB 
background. Throughout this Letter, we refer to these 
data samples as the ON and OFF samples, respectively. 

The Belle detector is a large- solid- angle magnetic spec- 
trometer described in detail elsewhere The main 
component relevant for this analysis is the electromag- 
netic calorimeter (ECL) made of 16.2 radiation lengths 
long CsI(Tl) crystals. The photon energy resolution is 
about 2% for the energy range relevant in this analysis. 

The strategy to extract the signal 6 — s- 57 spectrum is 
to collect all high-energy photons, vetoing those originat- 
ing from TT*' and rj decays to two photons. The contribu- 
tion from continuum e+e^ ~* QQ {q — u, d, s, c) events is 
subtracted using the OFF sample. The remaining back- 
grounds from BB events are subtracted using Monte- 
Carlo (MC) distributions scaled by data control samples. 

Photon candidates are selected from ECL clusters of 
5x5 crystals in the barrel region (—0.5 < cos9 < 0.84, 
where 9 is the polar angle with respect to the beam 
axis). They are required to have an energy E* larger 
than 1.5 GeV, and 95% of the energy has to be deposited 
in the central 3x3 crystal array. We require isolation 
cuts to veto photons from bremsstrahlung and interac- 
tion with matter. The center of the cluster has to be dis- 
placed from any other ECL cluster with E > 20 MeV by 
at least 30 cm at the surface of the calorimeter, and from 
any reconstructed track by 3 cm, or by 50 cm for tracks 
with a measured momentum above 1 GeV/c. Moreover, 
the angle between the photon and the highest energy lep- 
ton in the event has to be larger than 0.3 radians at the 
interaction point. We veto candidate photons from tt'^ 
and rj decays to two photons by combining them with 
any other photon. We reject the pair if the likelihood of 
being a 7r° or 77 is larger than 0.1 and 0.2, respectively. 
These likelihoods are determined from MC and are func- 
tions of the laboratory energy of the other photon, its 
polar angle 9 and the mass of the two-photon system. 

In order to reduce the contribution from continuum 
events, we use two Fisher discriminants. The first ex- 
ploits the spherical shape of BB events and is built using 
ten event-shape variables. These variables are calculated 
using either all tracks and showers in the event or ex- 
cluding the photon candidate. The event shape variables 
include Fox-Wolfram moments, thrust, and the angles 
of the thrust axis with respect to the beam and photon 
direction. The second discriminant exploits the topol- 
ogy of 6 — » 57 events and combines three energy flows 
around the photon axis. These energy flow variables are 
obtained using all particles, except for the photon candi- 
date, whose direction lies in the three regions defined by 
a* < 30°, 30° < a* < 140°, a* > 140°, where a* is the 
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FIG. 1: Photon energy spectra in the T(4S') frame. 

angle to the candidate photon. 

To optimize these selection criteria, we use a Monte 
Carlo simulation 0| containing large samples of BB, qq 
and signal weighted according to the luminosities of the 
ON and OFF samples. The signal MC is generated as 
a weighted sum oi B ^ K*j decays, where K* is any 
known spin-1 resonance with strangeness S = 1. The 
relative weights are obtained by fitting the total pho- 
ton spectrum to a theoretical model 0. The signal MC 
is normalized to the average measured branching frac- 
tion 1^. To improve the understanding of the photon 
energy spectrum at low energies, the selection criteria 
are optimized to maximize the sensitivity to the signal in 
the energy bin 1.8 GeV < E* < 1.9 GeV. 

After these selection criteria we observe 1.2 x 10^ pho- 
ton candidates in ON and 1.1 x 10^ in OFF data. The 
spectrum measured in OFF data is scaled by luminosity 
to the expected number of non-BB events in ON data 
and subtracted. To take into account the effect of the 
60 MeV (0.5%) energy difference, the measured OFF en- 
ergies are scaled by an empirical factor of 1.004 obtained 
from a MC study. The ON and scaled OFF spectra and 
their difference are shown in Fig. ^ 

We then subtract the backgrounds from B decays from 
the obtained spectrum. Five background categories are 
considered: (i) photons from tt^ — > 77, which account 
for more than half of the background in the 1.8-2.8 GeV 
range; (ii) photons from rj — > 77; (in) other real photons 
(mainly decays of w, 77', and J/ip, and bremsstrahlung); 
(iv) ECL clusters not due to single photons (mainly elec- 
trons interacting with matter, and n); (v) beam 
background. For each of these categories we take the 
predicted background from MC and scale it according 
to measured yields wherever possible. The inclusive 
B TT^X and B rjX spectra are measured in data 
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using pairs of photons with a well-balanced energy and 
applying the same ON— OFF subtraction procedure. The 
yields obtained are 5 to 15% larger than MC expecta- 
tions depending on the photon energy range. Since there 
is good agreement between MC and data for all features 
of the GEANT simulation for photons and electrons, we 
believe that the observed discrepancy between the mea- 
sured and simulated 7r° spectrum is due to the genera- 
tor Beam background is measured using a sample 
of randomly triggered events and added to the BB MC. 

For each selection criterion and each background cate- 
gory we determine the i?*-dependent selection efficiency 
in OFF-subtracted ON data and MC using appropriate 
control samples. We then scale the MC background sam- 
ple according to the ratio of these efhciencies. The effi- 
ciencies of the tt" and rj vetoes for non-Tr", non-r/ photons 
are measured in data using one photon from a well re- 
constructed 7r° applying the veto without using the other 
photon of the pair. The tt'' veto efficiency is measured 
using a sample of photons coming from measured 
decays. We use partially reconstructed Z)*+ D'^tt^, 
K~TT^Ti^ decays where the 7r° is replaced by the 
candidate photon in the reconstruction. The rj veto ef- 
ficiency for photons from vr'^'s and event-shape criteria 
efficiencies are measured using a tt*^ anti-veto sample. It 
is made of photons passing all selection criteria except 
the 7r° veto, which are combined with another photon in 
the event to give a Tr^-likelihood larger than 0.75. Other 
efficiencies are measured using the signal sample. 

The ratios of data and MC efficiencies versus E* are 
fitted using first or second order polynomials, which are 
used to scale the background MC. Most are found to 
be statistically compatible with unity. An exception is 
the efficiency of the requirement that 95% of the energy 
be deposited in the central nine cells of the 5x5 clus- 
ter, which is found to be poorly modelled by our MC 
for non-photon backgrounds. We estimate the efficiency 
for data using a sample of candidate photons in OFF- 
subtracted ON data after subtracting the known contri- 
bution from real photons. This increases the yield of 
background (iv) by 50%. The yield from the five back- 
ground categories, after having been properly scaled by 
the above described procedures, are subtracted from the 
OFF-subtracted spectrum. The result is shown in Fig.^ 

The spectrum contains 24100 ± 2140 ± 1260 events in 
the 1.8-2.8 GeV energy range, where the two errors are 
the statistical error of the OFF-subtracted ON data and 
of the BB background subtractions, and the systematic 
error related to the data/MC efficiency ratio fits used in 
the BB background scaling. We correct this spectrum 
for the signal selection efficiency function obtained from 
signal MC, applying the same data/MC correction fac- 
tors as for the generic photon background category (in). 
The average signal selection efficiency is 23%. 

The efficiency-corrected spectrum is shown in Figure|21 
The two error bars for each point show the statistical 



Source of systematic error 


xlO-* 


Raw branching fraction 


3.51 ±0.32 


Data/MC efficiency ratio fits 


± 0.208 


Choice of fitting functions 


± 0.048 


Number of BB-events 


+ 0.139 
- 0.160 


ON-OFF data subtraction 


±0.026 


Other BB photons 


± 0.054 


77 veto efficiency on 77 


± 0.008 


Signal MC 


± 0.089 


Photon detection efficiency 


± 0.072 


Energy leakage 


+ 0.035 
- 0.000 


Total error for partial B{b q-y) 


+ 0.282 
- 0.291 



TABLE I: Overview of systematic errors. 



and the total error, including the systematic error which 
is correlated among the points. As expected, the spec- 
trum above the cndpoint for decays of B mesons from 
the T(4S') at about 3 GeV, is consistent with zero. Inte- 
grating this spectrum from 1.8 to 2.8 GeV, we obtain a 
partial branching fraction of (3.51 ± 0.32 ± 0.29) x lO""*. 

The sources of systematic error are listed in Table U 
They are added in quadrature. The largest sources are 
the errors of the data/MC efficiency ratio fits (5.9% of 
the signal yield). For the error related to the choice of the 
polynomial functions in the data/MC efficiency ratio fits, 
we perform the same fit increasing the polynomial order 
by one. The number of BB events is determined from the 
number of hadronic events in ON and OFF data. The 
relative luminosities of the two samples are determined 
from radiative Bhabhas and e+e^ — > /^^M~ events. The 
errors on the OFF data subtraction arc estimated using 
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FIG. 2: Efficiency-corrected photon energy spectrum. The 
two error bars show the statistical and total errors. 



5 



the result of the fit to the spectrum above the endpoint. 
We integrate the resulting function in the 1.8-2.8 GeV 
range and obtain a yield of +40 ± 160. We add ±200 
(0.8%) to the systematic error. As we do not measure 
the yields of photons from sources other than tt'^'s and 
ry's in BB events, we vary the expected yields of these 
additional sources by ±20%. For the model dependence 
of signal selection efficiency we use an alternate signal 
MC that favors high-mass resonances decaying into high- 
multiplicity final states. Using this MC to correct for 
the efficiency changes the branching fraction by ±2.5%. 
The error on the photon detection efficiency in the ECL 
is measured to be 2.3% using radiative Bhabha events. 
This error also affects the estimation of photons from BB 
and contributes ±2.0% to the systematic error. Due to 
the low-energy tail in the photon energy measurement, 
some part of the spectrum may lie below the range of 
integration. We estimate this fraction to be smaller than 
1%. As this value is shape-dependent, we do not correct 
the measured branching fraction for it but instead add a 
(lao)% systematic error. 

In order to obtain the total 6 — > 57 branching fraction 
we apply corrections for the contribution from Cabibbo 
suppressed b dj decays and for the invisible part of 
the spectrum below 1.8 GeV. The ratio of the b ^ sj 
and b ^ d'j branching fractions is assumed to be Rd/s = 
(3.8 ± 0.6)% 0. The selection efficiency for b ^ d-f is 
found to be equal to the efficiency for b ^ sj within 10%, 
which we include in the systematic error. The fraction 
of the spectrum above 1.8 GeV is assumed to be Ri,s = 
0.952 1q:o29 from Gambino and Misiak As a cross- 
check we also use the value from Kagan and Neubert 
Ri.s = 0.958 1 ^:^29 : and Ri.s = 0.95±0.01 from Bigi and 
Uraltsev 0. We combine the errors on Rd/si Ri.s and 
the difference between the Ri,s values into the theoretical 
error. With these two corrections, we obtain 

B{b^ S7) = (3.55 ± 0.32 1 o fi - 0.07) x 10"^ 

for the total b ^ s"f branching fraction. This result is in 
good agreement with theoretical expectations and with 
previous experimental measurements [kI Il6j . 

We also measure the first two moments of the energy 
spectrum in the B rest frame. We extract the raw mo- 
ments from the distribution shown in Fig. |21in the range 
1.8 GeV < E* < 2.8 GeV and correct them for the ef- 
fect of the boost of the B meson in the T{AS) frame, for 
the energy resolution and for the 100 MeV binning. We 
do not correct the moments for the missing low-energy 
tail. We obtain the following moments for i?* > 1.8 GeV 
(corresponding to > 1.815 GeV in the B rest frame): 

{E^) = 2.292 ± 0.026 ±0.034 GeV 
(E^) - {E-yf = 0.0305 ± 0.0074 ± 0.0063 GeV^ 

where the errors are statistical and systematic. 



The systematic error contains the errors related to the 
moments corrections and the error sources already men- 
tioned for the branching fraction extraction. For the 
first moment the systematic error is dominated by the 
data/MC efficiency ratio fits (±0.9%) and the shape of 
the energy resolution (±1.0%). The error on the sec- 
ond moment is dominated by the data/MC efficiency ra- 
tio fits (±17%). These results agree within la with the 
only previous measurement, done by the CLEG collabo- 
ration However, it should be noted that the CLEO 
results are obtained for E* > 2.0 GeV. 

In conclusion, we have measured the branching frac- 
tion and photon energy spectrum of 6 — s- 57 in the en- 
ergy range 1.8 GeV < -E* < 2.8 GeV in a fully inclusive 
way. For the first time 95% or more of the spectrum is 
measured, allowing the theoretical uncertainties to be re- 
duced to a very low level. Using 140 fb~^ of data taken 
at the T(4S') and 15 fb^^ taken below the resonance, 
we obtain B {b ^ 57) = (3.55 ± 0.32 1 ° j° 1 ° x 10-^ 
where the errors are statistical, systematic and theoret- 
ical, respectively. This result is in good agreement with 
the latest theoretical calculations [3, ■ We have also 
measured the moments of the distribution and obtain 
(Ey) = 2.292 ± 0.026 ± 0.034 GeV and (E^) - {E^f = 
0.0305 ± 0.0074 ± 0.0063 GeV^ for E* > 1.8 GeV, where 
the errors are statistical and systematic. 
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